Accepted for publication in ApJ Letters 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



oo 



> 



CO 

o 

On 
On 

S3 

Oh: 

6 

CO 



X 



THE X-RAY LUMINOSITY FUNCTION OF NEARBY RICH AND POOR CLUSTERS 
OF GALAXIES: A COSMOLOGICAL PROBE 

Michael J. Ledlow 1 , Chris Loken 2 , Jack O. Burns 2 , Frazer N. Owen 3 , 

and Wolfgang Voges 4 

institute for Astrophysics, University of New Mexico, Albuquerque, NM 87131 
2 Office of Research and Dept. of Physics & Astronomy, University of Missouri 
Columbia, MO 65211 
3 National Radio Astronomy Observatory^], Socorro, NM 87801 
4 Max Planck-Institut fur Extraterrestrische Physik, Postfach 1603, D-85740, 
Garching bei Miinchen, Germany 
Accepted for publication in ApJ Letters 

ABSTRACT 

In this letter, we present a new determination of the local (z < 0.09) X-ray luminosity function 
(XLF) using a large, statistical sample of 294 Abell clusters and the ROSAT All-Sky-Survey. Despite 
the optical selection of this catalog, we find excellent agreement with other recent determinations of the 
local XLF. Given our large sample size, we have reduced errors by s=s a factor of two for Lx(o. 5-2.0 keV) > 
10 43 /i^q 2 ergs/ sec. We combine our data with previous work to produce the most tightly constrained local 
determination of the XLF over three orders of magnitude in Lx in order to explore possible constraints 
imposed by the shape of the XLF on cosmological models. A set of currently viable cosmologies is used 
to construct theoretical XLFs assuming LocM p and a erg — fin constraint (from Viana & Liddle 1996) 
based on the local X-ray temperature function. We fit these models to our observed XLF and verify that 
the simplest adiabatic, analytic scaling relation (e.g. Kaiser 1986) disagrees strongly with observations. 
If we assume that clusters can be described by the pre-heated, constant core entropy models of Evrard 
& Henry (1991) then the observed XLF is consistent only with 0.1 < fio < 0.4 if the energy per unit 
mass in galaxies is roughly equal to the gas energy (i.e., if j3 ~ 1). 

Subject headings: cosmology: observations - galaxies: clusters: general - X-rays: general 



1. INTRODUCTION AND BACKGROUND 

Much of the work on the luminosity distributions of 
rich clusters has been motivated by the results of Henry 
et al. (1992) who found evidence for statistically signifi- 
cant negative evolution in the XLF (i.e. fewer high Lx 
clusters at higher z) at z > 0.3 for £x(o.3-3.5fcevn > 
5 x 10 44 /i^" 2 ergs sec -1 from 67 clusters in the Einstein 
Extended Medium-Sensitivity Survey (EMSS). Recently, 
Vikhlinin et al. (1998) have confirmed the EMSS result at 
z > 0.3 for Lx(o.5— 2keV) > 3xl0 44 /i^" 2 ergs/ sec from a 160 
deg 2 survey from pointed ROSAT fields. They found a fac- 
tor of 3-4 decrease in the number density of these high Lx 
clusters as compared to a zero-evolution model. Several 
other studies have claimed no evolution in the XLF out to 
redshifts as high as z = 0.8 (Burke et al. 1997; Jones et al. 
1998; Rosati et al. 1998). However, none of these studies 
have sufficiently large search volumes to address evolution 
in the XLF at the highest X-ray luminosities and thus do 
not contradict the original EMSS result. 

Of prime importance in any evolutionary study is an ac- 
curate determination of the local XLF as a baseline to com- 
pare with the distant cluster XLF. Until recently, even the 
local XLF was quite poorly constrained due to low cluster 
numbers. The largest local samples compiled to date are 
the X-ray Brightest Abell Clusters (XBACS) (Ebeling et 
al. 1993, 1996) and the Brightest Cluster Sample (BCS) 
of Ebeling et al. (1997, 1998). The BCS includes 199 X- 



ray selected clusters down to w 5 x 10 42 ergs sec -1 in the 
0.1 — 2.4keV band out to z < 0.3. Consistent with most 
previous claims, no evidence was found for evolution in the 
XLF within z < 0.2 - 0.3 (Ebeling et al. 1998). 

We have examined a statistically complete sample of 294 
Abell rich clusters within z < 0.09 using the ROSAT All- 
Sky-Survey (RASS) over the energy band 0.5 — 2 keV as 
part of a multiwavelength study of nearby galaxy clusters. 
Unlike most other studies, our sample is purely optically- 
selected within the criteria for inclusion in AbelPs catalog. 
There is some overlap with both the BCS and XBACs 
sample, with the primary differences that we have used 
only Abell's northern catalog (Abell 1958), and our X- 
ray flux- limit is approximately a factor of 8 lower than the 
BCS sample. Our sample is larger than the BCS while our 
volume is nearly 30 times smaller. Given our large sample 
size, we have reduced statistical errors in the local XLF for 
Lx > 10 43 /i^" 2 ergs sec -1 by up to a factor of 2 compared 
to previous work. Combined with the poor cluster XLF 
of Burns et al. (1996) (BLL96), we examine the composite 
local XLF over more than 3 orders of magnitude in Lx in 
order to understand the cosmological constraints imposed 
by the tight power-law shape noted in BLL96. 

In section 2 we describe the sample, the derivation of the 
local XLF, and discuss the limitations imposed by our sam- 
ple selection. In section 3 we compare our new XLF with 
previous work. In section 4 we explore the consequences of 
the shape of the local XLF with regard to Press-Schechter 
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analytic predictions of the mass-function and possible con- 
straints on H,q and A. We list our conclusions in section 
5. We adopt Ho = 50 hso km sec -1 and qo = 0.5 when 
dealing with the observational data. 

2. THE SAMPLE AND DERIVATION OF THE XLF 

Our cluster sample is derived from Abell's Northern cat- 
alog, and includes all Abcll clusters in the range 0.016 < 
z < 0.09 with galactic absorption less than 0.1 magnitudes 
at R-band (logN H « 20.73). See Voges et al. (1999) and 
Lcdlow & Owen (1995) for more details on the sample se- 
lection. The total sample includes 294 Abcll clusters. All 
clusters have measured redshifts and we include all rich- 
ness classes in the sample. We calculate a survey area of 
14,155 deg 2 or 34% of the sky. Within our observed vol- 
ume we find the number density of clusters to be constant 
as a function of richness class and redshift suggesting that 
our sample is nearly complete and volume-limited within 
the limits of Abell's selection criteria. These findings are 
consistent with those of Briel & Henry (1993) and Mazure 
et al. (1996) with regards to the completeness of Abell's 
catalog over this redshift regime. 

The X-ray luminosity function was derived from images 
produced by the RASS as described in Voges et al. (1999). 
X-ray luminosities were calculated within a metric aper- 
ture of 0.75 h§Q Mpc in diameter over the energy band 
0.5-2 keV assuming a thermal spectrum with T = 5 keV. 
Corrections for missing flux were made according to the 
prescription of Briel & Henry (1993) (using (3 = 2/3) to 
produce a total Lx for each cluster over our ROSAThand. 
The primary effect of using a different (3 would be to shift 
the total luminosities to higher or lower values (a larger [5 
results in a smaller correction, thus lower total Lx), while 
not significantly changing the shape or amplitude of the 
XLF within the error bars. 

Voges et al. found a total detection rate of 83% for this 
sample of Abell clusters. For non-detections, we adopt the 
3cr upper- limits given in their Table 1. Because of varia- 
tions in exposure time (and slight variations in galactic 
absorption) across the sky with the RASS, each cluster 
has a different flux- limit, or maximum volume to which 
the cluster could have been detected. We follow the pre- 
scription of Avni & Bahcall (1980), and calculate the ob- 
served volume separately for each cluster. The volume is 
evaluated from z mi „=0.016 to the maximum redshift at 
which the cluster could have been detected with a 3a con- 
fidence. For clusters with only upper-limits to Lx we set 
Zmax equal to the redshift of the cluster. The XLF is then 
found by calculating dn(L)/dL as the sum over all clusters 
divided by the maximum search volumes of each cluster. 
Each binned data point is then found by dividing the above 
sum by the binwidth (ALx). For the entire sample, we 
find (y/V m ax) = 0.56±0.02. Error bars on the data points 
were calculated assuming Poisson statistics following the 
prescription of Rosati et al. (1998). 

3. THE X-RAY LUMINOSITY FUNCTION 

In Figure 1 we show the differential XLF for our low- 
redshift cluster sample. Also on this plot are the mea- 
surements of BLL96 derived from 49 poor clusters and the 
BCS sample of Ebeling et al. (1998). The steady decline 
in volume-density observed in our rich cluster sample for 



Lx < 10 43 /i^q 2 ergs /sec can be understood from the lim- 
itations of Abell's optical selection criteria. Because Lx 
varies considerably for a given optical richness (Voges et 
al. 1999), there are a significant number of optically poor 
clusters with Lx in the range of Richness Class clusters 
which are not in our sample. Thus, our sample is truly 
volume-limited only for clusters above this cutoff in Lx- 
Note, however, that for Lx > 10 43 ergs/ sec, our Abell 
cluster sample and the BCS sample are in excellent agree- 
ment. The BCS also extends to higher Lx because of the 
larger search volume (z < 0.3). Our XLF shown in Figure 
1 is also consistent with that of Edge et al. (1990) and 
Briel & Henry (1993). 



_ .001 r m 
"C-0001 s- 

w Z 

a io- 5 r 

CD 

| IO" 6 r 

"o 10" 7 =- 
Q. E 

5 = 
£ IO" 8 r 

1 10 i- 
10-'° i- 



m 



.01 .1 1 

L x (0.5-2.04 keV) (IO 44 " -2 ergs sec"') 



10 



Fig. 1. The X-ray luminosity function derived from our low-rcdshift 
{z < 0.09) Abell cluster sample (solid circles), the poor-cluster data 
points from Burns et al. (1996) (open squares), and the XLF from the 
Brightest-Cluster Sample (BCS) of Ebeling et al. (1998) (open circles). 

The local, differential XLF is remarkably well repre- 
sented by a power-law over more than three orders of mag- 
nitude in Lx- The high luminosity break in the XLF oc- 
curs at > 10 45 ft-^" 2 ergs/ sec, and can be seen when we 
include the highest luminosity point from the BCS sam- 
ple. Using the combined XLF of BLL96 and our new 
determination of the local rich-cluster XLF (for Lx > 
10 43 /i^ 2 ergs /sec), we find a power-law fit of the form 
4>(L) — KL^ where L44 is the X-ray luminosity in units 
of 10 44 ergs sec 1 and K is in units of 10 7 Mpc 3 . 
We find best-fit values of a = 1.83 ± 0.04 and K = 
2.35±Q 22- For completeness, we also fit a Schechter func- 
tion after including the highest-Lx point from the BCS. 
For a fit of the form: ^ = A exp(-L x /L* x ) L x a , we 

find A=(2.93±0.14)xl0- 7 (Mpc- 3 L^ 1 ), L x{0 5 _ 2kcV) = 

5.49±0.39 (10 44 ergs/ sec), and a = 1.77±0.01. These val- 
ues are consistent within the errors to the BCS, the RDCS 
XLF (Rosati et al. 1998) out to z=0.6, and the Southern 
SHARC survey (Burke et al. 1997) for 0.3 < z < 0.7. Note 
that these results do not conflict with the claimed nega- 
tive evolution in the XLF observed by Henry et al. (1992), 
and most recently by Vikhlinin et al. (1998) at the highest 
luminosities. 

As noted by BLL96, the remarkable power-law shape 
over such a large range in Lx suggests a continuity in that 
the bulk X-ray properties of poor clusters must not be fun- 
damentally different from richer systems. We explore the 
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consequences of this result in the next section. 

4. DERIVATION OF THE THEORETICAL XLF 

In order to assess the constraints our local XLF imposes 
on cosmological models, we compare it with various ana- 
lytic predictions. We proceed by using the Press-Schechter 
(PS) formalism (e.g., Press & Schechter 1974; Bond et 
al. 1991) to construct theoretical mass functions and then 
convert these to XLFs assuming a form for the X-ray mass- 
to-light ratio (c.f., Evrard & Henry 1991; hereafter EH91). 

We begin with the set of cosmological models whose pa- 
rameters are listed in Table 1. These models form a repre- 
sentative sample of current views as they include open and 
flat universes spanning a range in Cl a . For each model, the 
rms density fluctuation on 8/i~ 1 Mpc scales (<j$) was deter- 
mined from the as — Q relation of Viana & Liddle (1996) 
which, in turn, was fixed by the local number density of 
7 keV clusters. The Hubble constant was chosen to give 
an age for the Universe of roughly 12.5 Gyrs (consistent 
with globular cluster age determinations; e.g., Chaboyer et 
al. 1998). For each model we list the relative contributions 
of matter (fl ), baryonic matter (fib), and the cosmolog- 
ical constant (ft\) to the overall energy density. Power 
spectra for all the models were generated using the code 
described in Klypin & Holtzman (1997) and then PS mass 
functions (with 5 C = 1.3) were computed at z = 0. 

Our PS mass functions can be converted to XLFs by 
assuming a form for the mass-luminosity relation and cor- 
recting to our bandpass. We assume the bolometric X- 
ray luminosity is related to cluster mass as Lboi — cM p 
and will later fit for the parameters c and p. There ex- 
ist at least two theoretical predictions for the value of the 
exponent p. The self-similar model of Kaiser (1986), de- 
rived assuming a power-law initial perturbation spectrum 
and purely adiabatic gas physics, predicts p = 4/3 but 
it is well known that this fails to give the correct shape 
for the XLF (e.g., EH91; see also below). However, pre- 
heating of the ICM at an early epoch (possibly by galaxy 
formation) results in a different scaling relation and also 
resolves several discrepancies between theoretical and ob- 
servational results concerning evolution in the XLF (e.g., 
Evrard 1990, Navarro, Frenk & White 1995). For the case 
of a constant entropy core, EH91 derived a scaling which 
implies p = (10/3 — 3)/3/3 where j3 = /j,m p a 2 /kT is the 
usual ratio of dark matter to gas 'temperatures'. 

We correct our bolometric luminosities to the 0.5-2 keV 
bandpass by calculating temperatures and applying a cor- 
rection appropriate for a plasma with a metallicity of 
Z=0.3Zq. Specifically, the temperature corresponding to 
a given mass can be calculated from the analytic M-T re- 
lation derived from the virial theorem (e.g., Bryan & Nor- 
man 1998): kT = ( 1Ql f M0 Al /3 h 2 ' z keV where 

A c is the current density contrast within the cluster virial 
radius. The luminosity in our bandpass is then calculated 
by applying the usual bremmstrahlung correction factor as 
well as a multiplicative factor to account for the presence 
of metals (Bryan & Norman 1998; eqn. 21). 

Using the relation for Lboi and the bandpass correction, 
we converted our PS mass functions to differential lumi- 
nosity functions and made x-squared fits to a subset of 
the observational data. The observational points used in 
the fits are all four poor cluster points (BLL96), the five 



highest luminosity Abcll cluster points and the highest lu- 
minosity BCS point from Figure 1. We first set (3 = 1 in 
the M — T relation and fit for c and p. The fitted value 
for p is included in Table 1 and examples of two of the 
fits are shown in Figure 2. The dashed curve in Figure 
2a is the best fit when the exponent is kept fixed at the 
analytic prediction p — 4/3. Clearly, the shape of the XLF 
derived using this prediction is in gross disagreement with 
the observed function. Figure 2b also shows the impor- 
tance of the low and/or high- luminosity data points. If 
only our five Abell cluster data points are used (dotted 
line), the fitted value of p increases by at least 0.2 in all 
cases (from p = 3.18 to p = 3.88 in this case). We get vir- 
tually identical results if we redo our fits without the BCS 
point whereas dropping the poor cluster points results in 
slightly greater discrepancies. 
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Fig. 2. Fits to the observed XLF for two models. The solid lines are the 
best fits to a subset of the observational XLF points for models ACDM1 
(left panel) and ACDM5 (right panel). The dashed line in the left panel 
is the best-fit to the data assuming p— 4/3 as predicted by the self-similar 
analytic models. The dotted line in the right panel represents the best-fit 
to the data when the BCS and poor cluster points are ignored. 

If we invoke the constant entropy core model of EH91, 
then the exponent in the mass-luminosity relation is ac- 
tually a function of /3 (p = (10/3 — 3)/3/3). In this case, 
we fit for c and [3 and find the values listed in Table 1. 
Interestingly, only the models with 0.1 < fio < 0.4 are 
consistent with the expected value /3 ~ 1 . [A recent obser- 
vational analysis found /3 = 0.94 ± 0.08 (Lubin & Bahcall 
1993) which is in good agreement with numerical results 
(e.g., Eke, Navarro, & Frenk 1998).] Thus, if the constant 
entropy core model of EH91 applies, the present-day XLF 
observations suggest a low-density universe but cannot dis- 
tinguish between open and flat cases. 

5. CONCLUSIONS 

Starting from an optically-selected, statistical sample of 
Abell clusters, we have made a new determination of the 
local XLF to compare to previous work and more distant 
cluster samples. Our cluster sample is larger than all pre- 
vious studies, and is contained within a smaller volume. 
For this reason, we have reduced statistical uncertainties 
in the local XLF by nearly a factor of two for a limited 
range in Lx [Lx > 10 43 Ii^q ergs /sec). It is only for 
Lx < 10 43 ergs/ sec that incompleteness due to the opti- 
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Table 1 

cosmological model parameters 



Model 



Aec a 



9-,, 



CT8 



OCDM1 
OCDM2 
OCDM3 
OCDM4 
OCDM5 
OCDM6 
ACDM1 
ACDM2 
ACDM3 
ACDM4 
ACDM5 



12.5 
12.7 
12.2 
12.7 
12.3 
11.6 
13.1 
11.8 
12.4 
12.5 
12.2 



70 
65 
65 
60 
60 
60 
80 
80 
70 
65 
60 



0.1 
0.2 
0.3 
0.4 
0.5 
0.7 
0.2 
0.3 
0.4 
0.5 
0.7 



0.026 
0.030 
0.030 
0.035 
0.035 
0.035 
0.020 
0.020 
0.026 
0.030 
0.035 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.8 
0.7 
0.6 
0.5 
0.3 



1.467 
1.162 
1.004 
0.897 
0.817 
0.702 
1.478 
1.160 
0.901 
0.864 
0.719 



1.7B1-2 
2.20^: 9 , 
2.50 
2.69 
2.86 
3.22 
2.14 _ 

2.49 j:?* 

2.77 
2 

3.18 



0.661 .05 



0.88^ 



+ .24 
.20 
+ .25 
-.22 
,+ .32 
- .21 
+ .29 
-.29 
+ .18 



+ 27 
-.22 
+ .31 
-.21 
3+31 



27 



'17 

D -.10 

1 14 + 41 

± - 1 ^ -.18 

1.47 tit 
1.90 t f 

4.o tii 

0-S4t.lt 
1 17 +- 33 

± - 11 -.19 

1.62 t f T 

i on +- 75 
±au -.50 

3 5 + 2 4 
JJ -1.6 



a current age of universe in Gyrs 

b Hubble constant in units of km/s/Mpc 

cal selection of our sample is apparent. The observed in- 
completeness is not a failing in Abell's catalog, but rather 
results from the contribution of poor clusters and groups 
below Abell's richness limit. 

Combined with the poor-cluster XLF of BLL96, we have 
examined the local XLF over nearly three orders of magni- 
tude in L x ■ We find that the local XLF is remarkably well 
represented by a power-law over nearly this entire range in 
Lx- This is significant evidence that hierarchical forma- 
tion results in similar cluster properties over a large range 
in Lx and mass. Including the brightest Lx clusters from 
the BCS sample which fall above the break in the XLF at 
Lx > 10 45 ergs /sec, we also performed a Schechter- 
function fit which is in good agreement with other recent 
surveys to much higher redshift (z < 0.7), confirming a 
lack of significant evolution at these luminosities. 

We have used our new local XLF to derive a constraint 
on Qq. This would appear to contradict a common claim 
that the eg — f^o dengeneracy can be broken only by includ- 
ing the evolution with redshift (e.g. Bahcall & Fan 1998). 
In fact, PS mass functions for combinations of a$ and Qq 
that satisfy a o~g — Qq constraint differ in shape. Borgani et 
al. (1999) have recently used the shape of the local XLF in 



order to constrain ct§ — ^o an d the shape of the L-T rela- 
tion. Including clusters at higher redshift, they concluded 
that flo — 0.4+Q 2 for open models, and flo < 0.6 for flat 
models assuming no evolution in the L-T relation; both 
of which are consistent with our results. In this work, we 
have used the shape of the local XLF, the local number 
density of 7 keV clusters, and the PS formalism in order 
to constrain the cluster M-L relation; Lx oc M p . There is 
a clear trend for p to increase with fio (see also Mathiesen 
& Evrard 1998). None of the theoretical models are con- 
sistent with the analytic prediction p — 4/3 from Kaiser 
(1986). If we adopt the constant core-entropy model of 
EH91, and the additional constraint that /? w 1, the shape 
of the local XLF suggests that 0.1 < Slo < 0.4, with no 
constraint on A. 
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